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Abstract 

We define a CP sensitive asymmetry in the sfermion decays / — ► / Xj £ & , f Xj Q 9; 
based on triple product correlations between the momenta of the outgoing fermions. 
We study this asymmetry in the MSSM with complex parameters. We show that 
the asymmetry is sensitive to the phases of the parameters fi and M\. The leading 
contribution stems from the decay chain / — ► / x°j ~ > f Xi % ~~ * / Xi & & (/ Xi Q q)j 
for which we obtain analytic formulae for the amplitude squared. The asymmetry 
can go up to 3% for / — > / x± & A and up to 20% for / — ► / x? Q Q- We also estimate 
the rates necessary to measure the asymmetry. 



1 Introduction 



In the standard model (SM) the source of CP violation is given by the phase in the 
Kobayashi-Maskawa matrix [1]. However, it has been argued that this source is not 
enough to explain the observed baryon asymmetry of the universe (see for example [2]) 
and new sources of CP violation have to be introduced. In the minimal supersymmetric 
extension of the SM (MSSM), several supersymmetric (SUSY) breaking parameters and 
the higgsino mass parameter can be complex. 

The phases of the SUSY parameters are restricted by the experimental upper limits 
on the electric dipole moments (EDMs) [3] of electron, neutron and the 199 Hg and 205 T1 
atoms. The limiting bounds are \d e \ < 4.3 x 10~ 27 ecm [4], \d n \ < 6.3 x 10~ 26 ecm [5], 
\du g \ < 2.1 x 10~ 28 ecm [6] and \dxi\ < 1.3 x 10 _24 ecm [7], respectively. The general 
consensus is that one of the following three conditions has to be realized: i) The phases 
are severely suppressed [8, 9], ii) supersymmetric particles of the first two generations are 
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rather heavy, with masses of the order of a TeV [10], iii) there are strong cancellations 
between the different SUSY contributions to the EDMs [11]. At one-loop level, for the 
electron EDM these are the neutralino and chargino contributions, and for the neutron 
EDM in addition also the gluino exchange contributes. While the phase of fi is restricted, 
Iwl ~ ^/lO, there is no such restriction on the phase of M 1 [12]. 

In order to clarify the situation an unambiguous determination of the SUSY phases 
is necessary. In particular, for determining also the sign of the phases, measurements of 
CP sensitive observables are necessary. The SUSY phases give rise to CP-odd (T-odd) 
observables already at tree level [13, 14, 15, 16]. An important class of such observables 
involves triple product correlations [17, 18]. They allow us to define various CP asymme- 
tries which are sensitive to the different CP phases. These observables could be measured 
at future linear collider experiments [19] and would allow us to independently determine 
the values of the phases. 

In this paper we consider a T-odd correlation in the decays 

with £ — e, /i, t, and q denotes a quark. The T-odd correlation for the leptonic decay is 
defined as 

Oodd = P/ • (Pt x Pi), (2) 
and that for the hadronic decay as 

O q odd = P f ■ (p q x p ff ), (3) 

where p denotes the three-momentum of the corresponding fermion. We define the cor- 
responding T-odd asymmetries as 

A e, q = r(Qf d g d > 0) ~ Wtk < °) = JSign[Q^|A^|MLips 

T r(cft d > o) + r(o& < o) ~ / |A^fdLi P s ' 1 ) 

where Ai e,q is the matrix element for the decay (1). For the measurement of A^ or A^ 
it is necessary to be able to distinguish between the charges of £ + and £~ or q and q. In 
the case £ — e, /i,t this should be possible experimentally on an event by event basis at 
an e + e~ linear collider [19]. A\ may be measureable in the case of q — c, b, where flavour 
reconstruction is possible [20]. However, this will only be possible statistically for a given 
event sample. 

The leading contribution to the triple products (2) and (3) originates from the decay 

chain 

/ - fx] ^fxiZ^ fxl £ £ (fXi q ?), (5) 

which is shown schematically in Fig. 1. Essentially, the triple products (2) and (3) are 
correlations between the polarization and the Z boson polarization, which are encoded 
in the momentum vectors of the final leptons or quarks. The correlations would vanish, for 
example, if a scalar particle in place of the Z boson is exchanged. Final state interactions 
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Figure 1: Schematic picture of the subsequent two-body decays / — > fx}, x] ~~ * 
ZXi, Z — > U (qq) in the / rest frame. 

may also contribute to A^ 9 , however, they arise only at one-loop level. We expect that 
such contributions are smaller than 10%, because only weak corrections to the absorptive 
part of the x]-Z-Xi vertex have to be included. In the similar case of the decay xf —> 
W ± Xi, corrections smaller than 10% have been obtained [21]. Corrections of this order of 
magnitude have also been found in [22], where next-to- leading order effects on polarization 
observables within the SM have been studied. 

As will be shown, the tree-level contribution to the triple product correlations (2) 
and (3) are proportional to the imaginary part of the x^}-Z-Xi coupling squared and are 
sensitive to the phases of the neutralino mass parameters Mi and \i (see Eqs. (25)- (35)). 
O^dd is n °t sensitive to the trilinear scalar coupling parameter Af of the sfermion /. The 
reason is that the first decay in the chain (5), / — > fx}, is a two body decay of a scalar 
particle. In order to be sensitive to the phase of Af one would have to construct instead 
of (2), (3) a triple product correlation involving the transverse polarization of the fermion 
/. In principle this could be possible for / = r,t [23], but this case will not be considered 
here (for a similar case see [15]). The T-odd correlation (2) was proposed in [16] and the 
size of the asymmetry was calculated for the decay fi — > fi x\ ~ > Xi f 1 ^ + &~ ■> however, for 
a specific final state configuration only. In the present paper we extend the work of [16] 
by calculating the asymmetries (4) in the whole phase space. 

In Section 2 we give the definitions and the matrix element of the decay we are 
interested in. In Section 3 we perform the calculation of the T-odd asymmetry. In Section 
4 we present our numerical results. Section 5 contains a short summary and conclusion. 
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2 Definitions and Formalism 



2.1 Lagrangian and Couplings 

The parts of the interaction Lagrangian of the MSSM relevant for decay (5) are (in our 
notation and conventions we follow closely [24, 25]) 



where 



with 



£ /7x° = fkMj P L + 4jPR)x°j + h - c - > 3 = 1, • • • , 4 , k = 1, 2 , (6) 



4 = ^Lr4- b L = 9(KL) m BU> ( n = L > R ) ( 7 ) 



/lj = -^(tanOw^ + A^), 

f Rj = yfitwGwN^, (9) 



and 



Kj = (h u Rj )*=Y u N* 4 , 

t Lj = ^(tane^-2 + ^i), 

= -^-tanO^, (10) 



h d Lj = (h d Rj T = Y d N; 3 , 



V2 



f d Rj = ^-tanG^, (11) 

me r ] m„ 

Y td = — — l A ; y u = -= . (12) 

V2mvi/cos/9 \'2mw sin (3 

Here, Pl,r = 1/2(1 =F 75), g denotes the weak coupling constant, 0^ is the weak mixing 

angle, mw is the mass of the W boson and VJ kn is the scalar fermion mixing matrix defined 
below. m e and m u {m d ) is the mass of the corresponding lepton and up-type (down-type) 
quark, respectively. is the complex unitary 4x4 matrix which diagonalizes the neutral 
gaugino-higgsino mass matrix Y a p, N* a Y a pN£p = m^o5 ik , in the basis (B, W 3 , H®, H%) [24]. 
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The masses and couplings of the fk follow from the 
in the basis (/l, /r) reads 



£-m — (/L f- 



RJ 



{ e iLp f\M 2 ; I 

v JLR 



hermitian sfermion mass matrix which 

(13) 



fRR 



with 



ILL 

m 

JRR 



JUL 



M lf + ( J 4 _ 9/ sin2 ®w) cos 2/? ra| + m 2 f , 

M Rf + 9/ ®^ COS 2 ^ m l + m /' 

= axg[A f - pi* (cot fi) 21 ^], 



(14) 
(15) 

(16) 
(17) 



y^arg^-^cot/?)^], (17) 

where is the mass of the Z boson, q$ and l( L R is the charge and the isospin of the 
fermion, respectively. M L j, M R j,Af are the soft SUSY-breaking parameters of the fa 

system. The / mass eigenstates are (/i, / 2 ) = (/l, /r)^ with 

/ e iv f cos 9; sin 6 f ~ 



and 



cos 9 j = 



n r = 



JLR. 



The mass eigenvalues are 



J\Mj | 2 + (m| — Mj ) 2: 

V 7ifl 71 /LL 



\ — sin 0^ e l</3 / cos 0^ 

, sin (9? = = ■ 

f /l Mj 2 + (rrA - M 



/LL 71 



2 \2 



The remainin 



m A,2 = ^ (( M L + M L.) =f /P!^nJ) • 



(18) 



(19) 



(20) 



(5) are 
and 

respectively, where 



_ JLL JRR JLR / 

ig parts of the interaction Lagrangian of the MSSM relevant for the decay 
C zm ° = Z^(OfP L + Of P R )x) , h J = 1, • • • ,4 , 
C Zff -=Z»f lfl (L f P L + R f P R )f , 



(21) 
(22) 



0" L = 

lJ 4 cos W 



^—(N i4 N* 4 -N i3 N* 3 ) 
is Hw 



ij ~ ij 



(23) 

cosH,, " " ^ 
Note that our definition of 0'f ,R (Lf,Rf) differs from that given in [24] by a factor 
g/2cos<~V (9/ cos& w ). 



L f(Rf) = ~ /a (4l(r) - Qf sin 2 Q w ) . 

COS KJyy 
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2.2 Spin-Density Matrix Formalism 



For the calculation of the amplitude squared of the subsequent two-body decays (5) of 
the sfermion, we use the spin-density matrix formalism [26, 27]. The amplitude squared 
is given by 



\M\ 2 = |A(x?)| 2 \A(Z)\ 2 £ ^ { PDl ) XiX , ( PD2 )l\ ( PD3 ) 

A^A^A^A^ 



(25) 



with the propagators A(x°) = l/\p Xj — m^, + im Xj Y Xj \ and A(Z) = l/[p| — m| + 
irrizTz]- Here, p Xj ,m Xj , F x . (p z , m z , T z ) are the four-momenta, masses and widths of the 
decaying neutralino (Z boson), respectively The amplitude squared is composed of the 
unnormalized spin density matrices p^, po 2 and po A of the decay (5), which carry the 
helicity indices Aj, A- of the neutralinos and/or the helicity indices A&, X' k of the Z boson. 
Introducing a set of polarization basis 4-vectors s£ (a = 1,2,3) for the neutralino 
which fulfill the orthonormality relations s£ • s^. = — 5 ab and s£ • p Xj = 0, the density 
matrices can be expanded in terms of the Pauli matrices: 
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(p 

(Pd 3 ) 



D2)\ k ,X' k - 
A' fc A fe _ 



E p, Ev 



(26) 
(27) 
(28) 



The polarization vectors e^ k of the Z boson obey p z ^ k = and the completeness relation 
Y,\ k E u k * £ i k = — 5V + VzpVzvjirfz- The expansion coefficients of the density matrices 
(26)-(28) are 



D 



2 pa 



D^(\aL 2 + \blf)(p f .p Xj ), 



E^=^(|6j&| 2 -Kf)(P/'-^), 
4^ [2 Re(0'#0'« *)m Xl m X3 - (|0^| 2 + \0'{f\ 2 ){p xi ■ p Xj ) 
+ m'S\ 2 + \Olf\ 2 )(p X]P P Xl .+P X]a P XlP ) , 



(29) 
(30) 

(31) 
(32) 



yd 



i 8 ^MOfOf')^ fw-P** s xJ 

+*4w p^4>x,(l^l 2 + |o;'fi 2 ) 

p^m Xl Re{0'ifO'if*) , 

-2 <T (Lj + R 2 )(p f ■ p f ) + 2 (pf p*- + p£ #)(L* + it?) 
+i2(R 2 -L 2 ) p f , p fl/ , 
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(33) 



(34) 



with e 0123 = 1 and m Xl the mass of the lightest supersymmetric particle (LSP). The 
masses of the fermions / = e, /i, r, c, b are set to zero. In Eqs. (29) and (30) /' denotes the 
fermion steming from the first decay in Eq. (5). Inserting the density matrices (26)-(28) 
in Eq. (25), the amplitude squared is given by: 

\M\ 2 = 2 |A(x°)| 2 \A(Z)\ 2 {D 1 D 2pa + S« Di E^Df. (35) 



3 T-odd Asymmetry 

In the following we present in some detail the calculation of the T-odd asymmetry in 
Eq. (4) for the slepton decays £ — > £xj — > £Xi Z ~ * ^Xi f f ■ The replacements which 
must be made for the asymmetry in the case of q decays are obvious. 

In the rest frame of £, the coordinate system is defined such that the momentum 
four- vectors are given by 

p z = (.Ez,0,0, |pz|), p Xi = \p Xj \{E x ./\p x .\,wie 1 ,0,co8e 1 ), 

Pf — \Pf\(Ef/\Pf\, sin9 2 cos02, sin6> 2 sin0 2 , cos# 2 ) , (36) 

where 

IP/I - <VF__ l„ J mo fl.l ' W 



2 m~ t ; 2(E Z - \p z \ cos6 2 ) 



and 

± ( m l 3 ,+ m z ~ m l, ) I Pxj I cos ft ± ^ ^(m^ml.my -4|p Xj | 2 m| (1 -cos 2 #i) 
|Pzl = 2|p Xj | 2 (l-cos 2 6 1 ) + 2m 2 x . ' 

(38) 

with A(x, y, z) = x 2 + y 2 + z 2 — 2(xy + xz + yz). There are two solutions |pf | in the case 
|P X J < IPxj l' wnere |p°J = y / A(m 2 ^. , m|, m 2 J/2m^ is the neutralino momentum if the 
Z boson is produced at rest. The Z decay angle Q\ is constrained in that case and the 
maximal angle #™ ax is given by 

|Pv I mz A^(m 2 ., m 2 7 , m 2 ) 
sin gmax = ^ = ^ v x , < 1 . (39) 

|p x .| m z (m|-m 2 .) 

If |p° | > |p x J, the decay angle 6>i is not constrained and there is only the physical 
solution \p z \ left. 

The spin basis vectors of x*j i n the £ rest frame are chosen by 
4 =(o,^^l, 4. = fo, PxjXPZ , l, 4. = — (iPxI^xPx) , (40) 

Xj V |s 2 xs 3 |y Xj v Ipx^xp^I' Xj ™- v ' J j; 



™Xi 
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with p Xj = p Xj /\Pxj\- Together with p£./m Xj , the spin basis vectors form an orthonormal 
set. 

The Lorentz invariant phase space element for the decay chain £ — ^ £ x^ ^ £ Xi Z ^ 
i Xi f f, in the rest frame of £, can be written as 

dUps(m~,p e ,p xl ,Pf,p f ) = dUps(m~,p e ,p Xj )ds D2 

xJ2dUps(s D2 ,p Xl ,p^)ds D3 dUps(s D3 ,Pf,p f ) , (41) 
± 

where sc 2 = p Xj and sd 3 — p\. The Lorentz invariant phase space elements of the 
sequence of two body decays read 

dLips(mf,p,,p Xj ) = ^-y 2 [l - -2L J dQ , (42) 





Pz 


2 


\El\ 




COs6>] 




pIii 



dUps(s D2 ,p xl , Pz ) = , r-rrr dfti , (43) 



dLips( SD3 ,P/,P/) = si l riE l_^ lcos e 2) 2 d ^ > ( 44 ) 
where dfli = sin^ d6i d<pi. 

The partial £ decay width for the decay chain (5) is given by 

r(*~- l Xiff) = ^_J \M\ 2 dUps(mlp e ,p xl ,p f -,p f ). (45) 

We use the narrow width approximation for the propagators A(x°) and A(Z) : 

J \A(x°j)\ 2 ds D , 2 = — n r , / | A(Z)\ 2 dsD 3 = m ^r z - ^ ne approximation for the neutralino 

propagator is justified for (^-) 2 <C 1, which holds in our case with F Xj < 0(GeV). 
From Eqs. (4) and (35) we obtain for the asymmetry 

A t, q = I |A(x°)| 2 |A(Z)| 2 Sign[Q^ d ]S" Di S^gfdLips 

where in the derivation of this expression we have used / | A(x°)| 2 |A(Z)| 2 Sign[Oo' dd ] 
Di D 2pa D^ dLips = in the numerator and / |A(x°)| 2 |A(Z)| 2 E% 1 Ef^Df dLips = 
in the denominator. As can be seen from Eq. (46), the asymmetry A^ 9 is proportional to 
the spin correlation terms Ef^ Y^ a D . In the spin correlations only the term which 

contains the T-odd correlation e o '^ d , Eqs. (2) and (3), contributes to A^ q . The T-odd 
correlation O^d i s contained in the product of the first term of Eq. (33) and the last term 
of Eq. (34), which leads to 

Vd 2P *D? D -32 m xl lm(0';fO';f)(R} - L))e^ p XjP s« ja p f , p fv . (47) 



In the rest frame of £, (pe ■ s° ) = for a = 1,2, hence, = in Eq. (30) and only 
£|> 2 , defined in Eq. (33), contributes to the spin correlation terms in the total amplitude 
squared. Using the explicit representation of the fermion momentum vector, Eq. (36), and 
the neutralino spin vector, Eq. (40), the term with the e-tensor in (47) can be written as 

£ paiiv VxiP s Xj * Pfv Vjv = ~ m xj ■ (P/ x P/) = ~m Xj \p z \ |P/| sin^! sm6 2 sin0 2 , (48) 

where = p^/|p^|. Since < 9i, 62 < 7r the sign of the correlation p^ • (p/ x pj) is given 
by the sign of sin 02- Thus o dd > corresponds to an integration d02, while o dd < 
corresponds to an integration d02- We therefore integrate in Eq. (45) over the entire 
phase space except over the angle 2 - The T-odd asymmetry can then be written as 



A 



f 



' fir dr 

JO d<j> 2 


_ r2n dr " 

Jtt d<f>2_ 


d0 2 


' c-tt dr 

JO dfj> 2 


, r2vr dr " 
' Jt d<f>2_ 


d0 2 



(49) 



The dependence of A^ on the ik-^-Xj couplings on the Z-f-f couplings Lf,Rf 

and on the Z-x®-X°j couplings 0"j' R is given by 



which follows from Eqs. (46) and (47). Due to the first factor ~ J ~ J , the asymmetry 

\ a kj\ + \ b kj\ 

At will be strongly suppressed for the I ~ \b k A and maximally enhanced in the 



case of vanishing mixing in the slepton sector ; 3 ,„ ■ ) 3 in ±1. Due to the second factor 



L2 f- R2 f 

L'i+tii 



, is larger than with 



-4 W = f±f ffi^ffi -4 - 6-3 (4-5) x ^. (51) 



Note that the r.h.s of (47) and therefore the asymmetry in Eq. (4) vanish for m Xl 
0, which is related to the fact that it is possible to redefine the Weyl spinor xi e ta Xi 
in this limit. 



4 Numerical Results 



We present numerical results for the T-odd asymmetry A^ defined in Eq. (4). The 
values for A^ c may be obtained from Eq. (51). We analyze numerically the decay chain 
f — > tx® — > t~Xi Z — > TXi t £, £ = e, /i,t where Xi is the lightest supersymmetric 
particle (LSP). We assume that f x is the lightest sfermion and that the decays into a 
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real Xp 3 = 2, 3, and a real Z are kinematically possible. In the numerical study below 
we will treat the two cases f\ — > rx\ an d ^i — *■ T X3 seperately. The asymmetry A^ q in 
Eq. (4) could in principle also be studied in %° three-body decays if the two-body decays 
are kinematically forbidden [13], which will be treated elsewhere [28]. 

The relevant MSSM parameters are |/x|,</? M , \M\\, (pMn M 2 ,tan/3, \A T \, (pA T , m^, rrif 2 
and the Higgs mass parameter ttia- For all scenarios we keep tan/3 = 10, \A T \ = 
1000 GeV, (fiA T = 0, m fl = 300 GeV, m f2 = 800 GeV and use the GUT relation \M±\ = 
5/3 tan 2 QwM 2 in order to reduce the number of free parameters. We have checked that 
our results do not depend sensitively on this choice. We choose uia = 800 GeV to rule out 
decays of the neutralino into charginos and the charged Higgs bosons x°j -/> xfH^, i — 1,2, 
as well as decays into the heavy neutral higgs bosons 7^ Xi -^2,3 • For the calculation 
of the branching ratios BR(fx — > r x°) an d BR(x° — * Zxi), we take into account also 
the decays fi — > xJ^t, X°j ~^ H\Xa-> ^ ± Xa i n addition to f x — > x°r, x° — > ^Xi- -^1 i s the 
lightest neutral Higgs boson, which in general is not a CP eigenstate [29, 30, 31]. The 
decay f\ — > r£ £, £ — e, /i is kinematically forbidden due to our assumption that fi is the 
lightest sfermion. Other decay chains leading to the same final state are less important 
and will be neglected. 

4.1 Decay chain via X2 

First we consider A^ for the decay chain fi — > r x° ~~ r ^ X? — * T Xi^ ^ f° r £ = e, jj,, r. In 
Fig. 2a we show the contour lines for the branching ratio BR(ti ^ t Xi £ £) = BR(fi — > 
r x°) x BR(x° -> ^X?) x BR(Z -> £ I) in the M 2 -|/x| plane for = tt/2 and y? M = 0. 
Small values of the phase of \i are suggested by constraints on the EDMs for a typical 
SUSY scale of the order of a few 100 GeV [12]. For BRfa -> r ^ ^) we always sum 
over the lepton anti-lepton pairs which couple to the Z . The grey areas in Fig. 2 are 
kinematically forbidden since here m fl < m x o + m T (light grey) or m x o < m x o + mz (dark 
grey). We choose Mg > Mg since in this case the f\-T~x\ coupling \a[ 2 \ is larger, which 
implies a larger branching ratio BR(fi — > r X2) than for < Mg. (We use the usual 
notation = M Rf , M- L = M Lf , see Eqs. (14) and (15)). Mg > Mg is suggested in some 
scenarios with non-universal scalar mass parameters at the GUT scale [32]. Furthermore, 
in Eqs. (14) and (15) one could have M fLL < M fRR in extended models with additional 
D-terms [33]. In a large region of the parameter space we have BR(%2 —> ZXi) = I5 an d 
we take BR(Z — > £ £) = 0.1. The corresponding asymmetry A^ is shown in Fig. 2b. The 
dependence of A? on M 2 and is dominantly determined by lm(0"%0"2*)- 

In Fig. 3 we show the (fM 1 and ip^ dependence of BR(fx -» r xj f f ) and of A^ in 
the full range of the phases for \fi\ = 300 GeV and M 2 = 280 GeV. We display in Table 
1 the masses of Xi^ — 1> • • • A an d the total widths T X2 , T fl for various values of the 
phases. The value of A? depends stronger on </?Mi, which also determines the sign of A^, 
than on ip^. 

Based on our results on the asymmetry A? in fi — > r\\ — > Xi T ^ +1 ^~ an d the branch- 
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ing ratio we give a theoretical estimate of the number of produced fi's necessary to 
observe the T-odd asymmetry in Eq. (4). The relevant quantity to decide whether A^ 
is observable (at la), is given by ((A^) 2 x BR)" 1 , where BR stands for the branching 
ratio of the decay considered. The necessary number of produced fi's should then be 
> ((A^) 2 x BR)" 1 . As an example we take the point denoted by • in Fig. 3, with 
tp^ = tt/2 and = tt/2. For this point BR pa 2.5 x 10" 2 and \A e T \ pa 3 x 10" 2 which 
implies that {{A^) 2 x BR)" 1 pa 4.4 x 10 5 . For the decay fi — > &&x?r, on the other hand, 
BR pa 3.6 x 10" 2 and |^| pa 1.9 x 10" 1 , so that ((y^) 2 x BR)" 1 pa 7.7 x 10 2 . For 
comparison we consider a further example with smaller CP violating phases ip^ = and 
ipMi — — 0.37T (denoted by <g> in Fig. 3). Also in this case we obtain almost the same results 
for ((A e ^ b ) 2 x BR)" 1 . In these two examples, the asymmetry A^ q should be measurable 
at an e + e~ linear collider with y/s = 800 GeV and an integrated luminosity of 500 /6" 1 
for rn-fj = 300 GeV. It is clear that detailed Monte Carlo studies taking into account 
background and detector simulations are necessary to get a more precise prediction of the 
expected accuracy. However, this is beyond the scope of the present paper. For a Monte 
Carlo study on a T-odd observable in neutralino production and decay see [34]. 

Table 1: Masses of x®,i = 1, ... ,4 and widths F^ 2 , for various phase combinations 
of w and <p Ml , taking = 300 GeV and M 2 = 280 GeV, tan (3 = 10, A T = 1000 GeV, 
m fl = 300 GeV, m f2 = 800 GeV for M s > M L . 









X2i m X-ii m Xi 


[GeV] 


Tx 2 [MeV] 


r fl [MeV] 








135, 


234, 


306, 


358 


4.06 


527 





7T 

2 


137, 


233, 


308, 


357 


1.79 


550 





71 


138, 


231, 


309, 


356 


0.09 


573 


7T 

2 





137, 


239, 


307, 


353 


5.43 


487 


7T 

2 


7T 

2 


138, 


238, 


309, 


352 


2.89 


511 


7T 

2 


TT 


137, 


237, 


311, 


351 


1.49 


529 


7T 





138, 


245, 


309, 


347 


7.25 


448 


7T 


7T 

2 


137, 


244, 


311, 


346 


5.78 


466 


7T 


7T 


136, 


243, 


313, 


345 


4.32 


484 
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Figure 2: Contour lines of the branching ratio for fi — > ^\rit and asymmetry defined in 
Eq. (4) in the M 2 -yU plane for (p Ml = n/2 and (p^ = 0, taking tan/3 = 10, A T = 1000 GeV, 
m fl = 300 GeV, m f2 = 800 GeV for Mg > Mg. The grey areas are kinematically 
forbidden since here m fl < mo + m T (light grey) or m o < m a + m z (dark grey). 




Fig. 3a 





Figure 3: Contour lines of the branching ratio for f\ — > Xi T ^ an d asymmetry defined 
in Eq. (4) in the ^ M r^ plane for = 300 GeV and M 2 = 280 GeV, taking tan/3 = 10, 
A T = 1000 GeV, m fl = 300 GeV, m f2 = 800 GeV for M s > M L . The points denoted by • 
and <S>, respectively, are for the theoretical estimate of the necessary number of produced 
fi's (see text). 
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4.2 Decay chain via X3 

Next we discuss the decay chain fi — > r Xs — > i~Z Xi — > t"Xi^ The two decays fi — > r x° 
and Ti — > t ^3 can be distinguished by measuring the r energy in the fi rest frame. In 
Fig. 4a we show the contour lines for the branching ratio BR(ti — > r x° ^ ■O = BR(f 1 — > 
r x°) x BR(x° -> Zx?) x BR(Z -> £ I) in the M 2 -|//| plane for </? Ml = tt/2 and y? M = 0. 
The area A (B) is kinematically forbidden since m^o < m^o + mz (m^ < m^o + m T ). 
The grey area is excluded since m x ± < 104 GeV. We choose < since the fx-r-xl 
coupling \a[ 3 \ is larger, which implies a larger branching ratio BR(fx — > r xi) than for 
M% > Mi. The total branching ratio is smaller than for the previous decay chain since 
BR(fi — > t xl) < -75(0.05) in the upper (lower) part of Fig. 4a. 

The corresponding asymmetry is shown in Fig. 4b. The asymmetry A^ vanishes 
on contours where either \a[ 3 \ = \b[ 3 \ or lm(0'( 3 01 3 ) = 0, see Eq. (50). On the one hand, 
along the contour line in the lower part of Fig. 4b we have | a[ 3 \ = \bl 3 \. On the other 
hand, along the contour line in the upper part of Fig. 4b we have \m{0'[ 3 0'{ 3 ) = 0. 
Furthermore, between the upper and the lower part of Fig. 4b (area A), there is a further 
sign change of lm(0" 3 0" 3 *). Concerning the first factor in (50), we remark that it 
increases for increasing M 2 and decreasing \fj,\. This behaviour can be understood by 
observing that for \p\jM2 — > the gaugino component of x 3 gets enhanced, resulting in 
l&Ial/Ksl - 0. 

In Fig. 5 we show the dependence of BR(ri — > r Xi £ an d of A^ on the phases 
ifM 1 and (pp, fixing = 150 GeV and M 2 = 450 GeV. For these parameters we display 
in Table 2 the masses of = 1, ... ,4 and the total widths r^ 3 , for various phase 
combinations. Note that maximal CP violating phases ip^, p>u x = ±vr/2 do not necessarily 
lead to the highest value ofA^ due to the complex interplay of the phases in lm(Oi 3 0'{ 3 ). 
The value of A^ depends stronger on tpMi, which also determines the sign of A^, than 
on (pp. Comparing Fig. 3b and Fig. 5b, one can see that both figures have in common 
the strong <pm x dependence, where in a good approximation the sign of A^ is sgn(A^) ~ 
sgn^Mi) in Fig. 3b and sgn(A^) ~ —sgn(ipMi) in Fig. 5b. This difference can be traced 
back to the different behaviour of \m.(0'{%0"2*) and lm(0'{ 3 0'{ 3 *). Moreover, in Fig. 5b 
two points of level crossing appear at approximately ip Ml ps ±0.957r, ip^ pa ±0.7n. 
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Table 2: Masses of x®,i = 1, ... ,4 and widths T% 3 , for various phase combinations 
of (p,, and cp Ml , taking = 150 GeV and M 2 = 450 GeV, tan/3 = 10, A T = 1000 GeV, 
m fl = 300 GeV, m f2 = 800 GeV for M & <M L . 
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Figure 4: Contour lines of the branching ratio for fx — > X? r ^ an d asymmetry A l T defined 
in Eq. (4) in the M 2 -/i plane for tpM 1 = n/2 and ip^ = 0, tan/3 = 10, A T = 1000 GeV, 



rrii 



300 GeV, m* 



800 GeV for < M^. The area A (B) is kinematically forbidden 



since m^o < m^o+mz [m^ < m^,o+m T ) 



The grey area is excluded since m x ± < 104 GeV. 
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Figure 5: Contour lines of the branching ratio for fi — ► Xit££ and asymmetry A l T defined 
in Eq. (4) in the (p M r¥n plane for \n\ = 150 GeV and M 2 = 450 GeV, taking tan/5 = 10, 
A T = 1000 GeV, m fl = 300 GeV, m f . 2 = 800 GeV for Mb < Mr. 



5 Summary and Conclusion 

We have considered a T-odd correlation and the corresponding asymmetry in the sequen- 
tial decay / — > /' y® ~^ f Xi Z — > /' Xi //■ The analytical expressions have been given 
in the density matrix formalism. The contribution to the T-odd correlation is induced by 
possible CP violating phases in the neutralino sector. 

In a numerical study of the decay f\ — > tx\ ff we have shown that the T-odd 
asymmetry considered can be of the order of a few percent for leptonic final states. The 
number of produced f's necessary to observe is at least of the order 10 5 , which may 
be accessible at future collider experiments. For a semi-leptonic final state like Xi T bb 
the T-odd asymmetry is larger by a factor 6.3. If the T-odd asymmetry A\ of such a 
semi-leptonic final state could be measured with similar accuracy the number of produced 
f's necessary to observe A\ is of the order 10 3 . 
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